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' Abstract 
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I Choosing the three phenomenological models of the dynamical cosmological term 

A, viz., A ~ (a/a)^, A ~ d/a and A ~ p where a is the cosmic scale factor, it 
has been shown by the method of numerical analysis for the considered non-linear 
differential equations that the three models are equivalent for the flat Universe k = 
and for arbitrary non-linear equation of state. The evolution plots for dynamical 
cosmological term A vs. time t and also the cosmic scale factor a vs. t are drawn here 
for k = 0, +1. A qualitative analysis has been made from the plots which supports 
the idea of inflation and hence expanding Universe. 

Key words: general relativity, cosmological parameter, inflationary cosmology. 
BAGS: 04.20.-q, 04.20.Jb, 98.80.Jk 



Corresponding author (E-mail: saibal@iucaa.ernet.in). 
Preprint submitted to Elsevier 20 June 2008 



1 Introduction 



The observations on Supernovae type la [Tll^ lead tlie scientists to the concept 
of accelerating Universe. The unknown type of energy responsible for this kind 
of acceleration is known as dark energy. Now, there are various types of models 
related to this dark energy [3!H] and for those models the expansion rate is 
also different. The so called cosmological constant may be one of them. For 
the explanation of accelerating Universe, A of dynamical character is preferred 
rather than a constant one. This character can naturally follow from non-linear 
effects of gravitational and scalar fields, involved in the physical model. In the 
paper of Ray, Mukhopadhyay and Meng [5], it was shown that among the 
dynamical models of A, the three types A ~ (a/a)^, A ~ a/a and A ~ p 
are equivalent, where a is the cosmic scale factor of the Robertson- Walker 
metric and p is the matter-energy density. They also analytically established 
a relationship between the parameters a, (3 and 7 of the respective models. In 
the present work emphasis has been given to show the equivalence of the same 
three models of A by using the method of numerical analysis of underlying 
non-linear Einstein equation for the flat Universe {k = 0) and for equation of 
state with arbitrary non-linear dependence on matter-energy density. Another 
aspect of this work is to make an attempt for visualizing the incidents that 
occurred during the very early stage of the Universe, especially the feature of 
inflation P71l8] . 



2 The Field Equations and General Results 



The Einstein field equations are given by 
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where A= A{t) is the so called cosmological constant. Here c, the velocity of 
light, is assumed to be unity in relativistic units. 



Now, let us consider the Friedmann-Lemaitre-Robertson- Walker metric 



ds^ = -dt^ + a{tf 



dr"^ 
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+ r^(dr + sin' 



where the curvature constant A; = — 1, 0, +1 for open, fiat and close models of 
the Universe respectively. 



2 



For the above spherically symmetric metric (2), the non-linear Einstein field 
equations reduce to the following two equations, respectively the Friedmann 
equation and the Raychaudhuri equation as 
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where the fiuid pressure p and energy density p are related in the form 
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P = (5) 
u! being the equation of state parameter of the polytropic equation of state. 

The energy conservation law can be written as 
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Differentiating equation (3) with respect to time t we get 
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With the help of equation (6) equation (7) reduces to 
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Now, from equation (4) with the help of equations (8) and (5), we get 
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2.1 A ~ {a/af 



Let us now consider the following dynamical model of A which is 
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Employing this equation (10) in equation (3) one immediately obtains 
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Now taking n — 1 and using equation (11), equation (9) reduces to the form 



a = -(3^+1) 
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2.2 A~d/a 



As a second dynamical model we now start with 



By the use of equation (13) in equation (3), we get the value for p as 
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For n = 1, by the use of equation (14), the equation (9) finally gives the value 
of d as 



a — 



[1 + Su) 



[(^ + l)/?-2] [ 



a k 
— + - 
a a 



(15) 



2.3 A~p 



Let us now take the third form of phenomenological A as 



A = SnG-fp 



(16) 



where the constraint is such that 7 > 0. Using the value of A given in equation 
(16), we obtain from equation (3) the value of p as 
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Using equation (16) and (17), we get from equation (9) the value of a for n 



as 
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In this connection it is to be noted here that Ray, Mukhopadhyay and Meng [5] 
have shown that the three forms A = 3a(d/a)^, A = (3{a/a) and A = SnG^p, 
as expressed in equations (10), (13) and (16), are equivalent for = 0. They 
found through analytical method that the parameters involved in the three 
dynamical relations are connected by 



3(/3w + /3-2) 1 + 7- ^""^ 

This means that it is possible to find out the identical physical features of 
others if any of those three phenomenological A relations is known. 

In this connection it is, however, interesting to note here that for linear de- 
pendence on energy density in equation of state (5), corresponding to n = 1, 
the equation (15) of Case-2. 2 and the equation (18) of Case-2.3 are equivalent 
for k = +1,0,-1 when the relation 



(1 + 3u) 
[{cu + 1)P - 2] 

holds good. 
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3 Graphical Presentation of the Results 



Now using the method of numerical analysis, let us try to study the variation 
of the cosmological parameter A and the scale factor a with time which are 
shown in the following plots. 

From a close observation of the graphical plots, viz., figures 1(a), 2(a) and 3(a), 
one can find out equivalence of the three models with respect to time variation 
of A while figures 1(b), 2(b) and 3(b) exhibit equivalence of the same three 
models with respect to time variation of a for A; = as obtained in analytical 
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method by Ray, Mukhopadhyay and Meng |5]. It can also be observed that 
figures 2(c) and 3(c) show equivalence of A ~ a/a and A ~ p models with 
respect to variation of A with time for k = +1. However, the behaviour of A 
with respect to time for A ~ (a/a)^ model is quite different. For A ~ (d/a)^ 
model we observe an abrupt increase of the cosmological parameter within 
a very short period of time and then a comparatively slower decrease of it. 
This abrupt rise of A may be interpreted as the driving force behind inflation 
because time variation of the scale factor a for the same values of k, u, a, 
P and 7 show a very sharp increase as depicted in flgures 1(d), 2(d) and 
3(d). Moreover, the sudden jump in the value of A exhibited in flgure 1(c) for 
k = 1 shows a clear indication of the role of the dark energy candidate A as 
repulsive pressure in connection to inflationary phase of the Universe Pl7f8] 
and also may be interpreted as a numerical manifestation of the idea that 
dark energy is responsible for making the Universe flat during inflation. This 
is because immediately after attaining a peak value, A has dropped down. So, 
it is quite natural to think of that a huge amount of dark energy was used 
up for triggering the exponential cosmic expansion as well as for removing 
the curvature of the Universe. Interestingly, it is to be noted here that this 
type of qualitative variation of A with time during inflation was unavailable 
in the analytical method by Ray, Mukhopadhyay and Meng [S] . However, due 
to lack of advanced computing facility, in the present investigation we could 
not quantify the time duration of inflation which has been suggested in the 
literature as 10~^^ to 10"^^ sec [9]. 

The present work shows that although phenomenological models do not orig- 
inate from any quantum fleld theory, yet at least in some cases they can 
successfully reflect the present cosmological picture. This is another positive 
side of this investigation. In particular taking the values of Qmo = 0.330 ± 
0.035. p!|2|3P] . the ranges of the values of model parameters ao,/5o,7o, corre- 
sponding to the observed accelerated expansion, are obtained as 0.635 > ao > 
0.705, 3.417 >Po> 4.674 and 1.739 > 7o > 2.389. 

Using the above-mentioned values of ^Imo, the ranges of the present values of 
the cosmological parameter Aq are obtained as 1 x 10^^^ s^^ - 2 x 10^^^ s~^, 
which agree with the results of Carmeli [10] and Carmeli and Kuzmenko |llj . 
where they obtain the value 1.934 x 10~^^ s~^. 



4 Conclusions 

In the present investigation, instead of flnding out exact solution of the or- 
dinary non-linear differential equation, the method of numerical analysis has 
been adopted for the three phenomenological models of A, viz., A ~ (d/a)^, 
K ^ a/a and A ~ p. However, the main idea of the article i.e. the study of 
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numerical solutions for cosmological problems in absence of particular or gen- 
eral analytic solutions is not a new idea. Interesting literature in this aspect 
are available, where numerical cosmology with the Cactus code [12] or testing 
the Cactus code on exact solutions of the Einstein field equations [13] have 
been performed. 

The dynamical nature of A and its time variation of the present investigation 
can find physical origin in non-linear effects of gravity (as it was the case for 
nonlinear E? term induced by vacuum polarization in one of the first infla- 
tionary models [2]) or of scalar fields, involved in the model (as it took place 
in the of self consistent inflation [T^16|17] ). The interesting features of this 
investigation can be put in the following way: 

(1) The time variation of A and a for /c = support the work of Ray, 
Mukhopadhyay and Meng [5] so far as the equivalence of the three chosen 
phenomenological models are concerned; 

(2) It has been possible to establish that A ^ d/a and A ^ p models are also 
equivalent for /c = 1; 

(3) Finally, a bonus obtained from the present numerical work is the qualita- 
tive visualization of inflationary scenario of the Universe through the variation 
of A. 
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Fig. 1. Case-2.1 related to model A = 3a(d/a)^ shows variation of cosmological 
parameter for A; = 0. 
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Fig. 2. Case-2.1 related to model A = 3Q(d/a)^ shows variation of scale factor for 
A; = 0. 
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Fig. 3. Case-2.1 related to model A = 3Q(d/a)^ shows variation of cosmological 
parameter for k = +1. 
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Fig. 4. Case-2.1 related to model A = 80(0/0)^ shows variation of scale factor for 
k = +1. 




Fig. 5. Case-2.2 related to model A = P{d/a) shows variation of cosmological pa- 
rameter for k = 0. 
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Fig. 6. Case-2.2 related to model A = (3{d/a) shows variation of scale factor for 
k = 0. 
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Fig. 7. Case-2.2 related to model A = f3{a/a) shows variation of cosmological pa- 
rameter for k = +1. 
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Fig. 8. Case-2.2 related to model A = (3{a/a) shows variation of scale factor for 
k = +1. 
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Fig. 9. Case-2.3 related to model A = SirGjp shows variation of cosmological pa- 
rameter for k = 0. 
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Fig. 10. Case-2.3 related to model A = SttG'jp shows variation of scale factor for 
k = 0. 




Fig. 11. Case-2.3 related to model A = SttG'jp shows variation of cosmological 
parameter for k = +1. 



Fig. 12. Case-2.3 related to model A = SttG'jp shows variation of scale factor for 
k = +1. 
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